POLYNOMIAL MODELS FOR INTEREST RATES AND STOCHASTIC 

VOLATILITY 



SI CHENG AND MICHAEL TEHRANCHI 

Abstract. This essay discusses a class of arbitrage-free interest rate models that has tractable 
bond price. The resulting bond price depends on some suitable factor process in a polynomial way. 
We also extend this class of polynomial model to the stochastic volatility case and they can produce 
tractable price for power options. The main result of this essay is that we have established the 
necessary and the sufficient conditions for these models to be arbitrage free. 
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1. INTRODUCTION AND NOTATIONS 

In this document, we denote r< as the spot interest rate, Pt(T) as the time t price of Zero Coupon 
Bond (ZCB) with maturity date T > t, {St}t>o as the price process of the underlying asset, Ct(T, K) 
as the time t price of call option with maturity T and strike K and Pt (T, 9) as the time t price of 
power options that pays Sj, (0 < 6 < 1) at maturity date T > t. We also assume in this article that 
the filtration & is generated by some multidimensional Brownian motion W = (W 1 , . . . , W m ). For 
no arbitrage reason, we assume that under the EMM measure, the stock price process St follows: 

dS t = S t a t .dW t 

Where the multi-dimensional process at denotes the volatility. We also denote Z t as the scalor 
factor process that is the solution of some suitable SDE and we are going to construct the models 
based on this factor Z t . 

Remark 1.1. The fact that the factor process Zt is a scalor process will simplify the notation and 
calculation quite a lot. In general, when Z t is multi-dimensional, we still have similar results but 
the notation and calculation will be much complicated. 

There has been various famous class of models that produce tractable ZCB prices Pt{T) and 
power option prices P t (t,9). Most of them takes the form of exponential polynomial type. For 
example, in the case of interest rate case, we can model the ZCB price process as: 

P t (T) = e~fo~ trt( - x ^ dx 

The above formula encodes the price of ZCB Pt(T) by process rt(x) which is known as the forward 
rate and satisfy the boundary condition 

P T (T) = 1 VT 

The famous CIR [2] and Vasicek [8] model will produce the forward rate in the form 

r t (x) = A{x)Z t + B(x) 

where A, B are some deterministic functions. In the above examples forward rates depends on 
factor process Z t linearly, and for quadratic type forward rates, we can see examples in Leippold 
and Wu [5]. Surprisingly, under some conditions, these are the only examples that produce no 
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arbitrage models. As Filipovic [3] points out that the maximum degree that the forward rate can 
have is 2, i.e. quadratic. 

In the case of stock market, instead of encoding the call option price process Ct(T, K), which is 
the natural choice, we choose to encode the power option price process Pt(T, 9). The reason behind 
this is that sometimes we may fail to have a closed form solution to the call price Ct(T, K) but we 
may have closed form power option prices Pt(T, 9) [1]. 

Remark 1.2. Notice that knowing the process {P t (T,8)} for all T > 0,0 < 9 < 1 is equivalent to 
specify the conditional distribution of St given J^t. Therefore it is equivalent to specify the process 
{C t (T,K)} since under EMM 

C t (T,K)=E[(S T -K) + \^ t }. 

More precisely, we have the following identities: 

S e = 9(1 -9) (S A K)K e - 2 dK for all S > 0, 0 < 9 < 1 

J o 

The above idendity shows that a power option with payout S T can be replicted by a portfolio of 
options with payouts St f\K = St — (St — K) + , and hence by portfolio of call options and the stock 
itself. Similarly: 

-i r+co Tsl—0—ix 

SAK=— S e+ix - — — dx for allO <9<1 

27T./-OC (x - i9)(x + i(l - 9)) 

and hence options with payoff St A K can be replicated by generalised power options with payouts 
S T cos(x log St) and S T sin(x log St) 

Here we model the power option prices processes Pt(T,9) as: 

p t (T,0) = Sfe-* 9 ( 1 - 0 )£~ tft ( x > 0 ) dx 

Therefore we encode the process Pt(T, 9) by the process ft(x, 9) which we called forward variance. 
Notice that the above re-parametrisation satisfies Pr(T, 9) = S T for all T, 9. Similar to the interest 
rate case, we have examples that the forward variance ft(x,9) is a linear or quadratic function on 
factor process Z t . 



Example 1 (Heston's Model. [4]). The model set up is 

dS t = S tV %dW? 
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dZ t = {Z- Z t )dt + y/ZtdWl 

Where W s and W 1 are two Brownian motions with correlation p and Z is some constant. As a 
result, the forward variance ft(x,9) takes the form that: 

ft(x,6) = A(x,9)Z t + B(x,9) 

for some deterministic functions A,B. Hence a linear example. 

Example 2 (Stein-Stein Model. [7]). With the same notation as the last example, the model set 
up is 

dS t = S t Z t dWf 

dZ t = -idW} + (a + (3Z t )dt 

with some constant a,/3,j. The forward variance is given by 

ft(x, 9) = Q 2 {x, 9)Z 2 t + Q 1 (x, 9)Z t + Q 0 (x, 9) 

For some deterministic functions Q%. Hence a quadratic example. 

Similar as interest rate case, Michael and I have shown that under some constraints, the maximum 
degree that forward variance can have is up to 2. i.e. quadratic (See appendix A). Therefore together 
with a time dependent version of Black-Scholes Model (corresponds to the case ft(x, 9) = A(x, 9). 
degree 0 example), we have seen all the no arbitrage models. 

So far, all the tractable models introduced above have a common form: the price process Pt(T) 
or Pt(T,0) is of the form of exponential polynomial. In the rest of this paper, we are going to 
investigate a new class of tractable models that takes the form of polynomial instead of exponential 
polynomial. 

The outline of the paper is as follows. In section 2, we introduce the polynomial interest rate 
model and provide a sufficient and a necessary condition for these model to admits no arbitrage. 
In section 3, we do a similar treatment to the stock market case. And finally in section 4, we give 
some further analysis and provide some examples. 
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2. POLYNOMIAL INTEREST RATE MODELS 

2.1. Model Setup. In this section, we are going to discuss a class of interest rate models that 
takes the form: 



P t (T) = J29i(T-Wi(Zt 



i=0 

Where gi,<fii are some deterministic functions and the scaler factor process Z t follows: 

dZ t = a(Z t )dt + b{Z t )dW t 

Remark 2.1. Notice that if we set 4>i(z) = z l , the resulting Pt(T) will be in polynomial form. 

The boundary condition Pt(T) = 1 can be made to hold true by setting go(0) = 1, <?i(0) = 
0 i = 1, 2, . . . , n and 4>o{z) = 1. Here follows, we will always use this type of boundary condition. 

Remark 2.2. We are safe to consider only boundary conditions of the type above. If there exists 
some functions <fri(z) that is independent of z, then we can absorb the constant into the coefficient 
function gi(x) and assume (fti(z) = 1. If however, all 4>i(z) depends on z, we then have: 



J>(o)<M20 = i 



i=0 



9o(z) 

Hence our polynomial model reads: 

n n 

P t {T) = Y,9i{T - t)<t>i{Z t ) = g 0 (T - t)<h{Z t ) + J>*(T - t)<f>i(Z t ) 
i=0 i=l 

9o(T-t) ^ / g 0 (T-t) \ 

= + h' MT - " V - -^«r) MZt) 

We then can re-define: 

9o(T-t) 



90 (0) 

9o(T-t) 



G 0 (T-t) := 

Gi(T-t):= 9i (T-t)[l- 

V 90 (0) 

This re-parametrization will satisfy the above type of boundary condition. 

We also assume that the spot interest rate process r t also depends on the factor process Z t . i.e. 
r t = R(Z t ) for some deterministic function R. Hence the discounted ZCB price process Pt(T) can 
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be expressed as: 



P t (T) = exp (-J r s ds\ P t (T) 

= exp (- jf R(Z s )ds S j f>(T - t)<f>i(Z t 



i=0 

2.2. No Arbitrage Condition. For consistency, we need: 

(A) : The process P t (T) is local martingale under some EMM Q VT 
For no arbitrage, we have: 

Theorem 2.1. Condition (A) holds if and only if 

i=0 i=0 ^ ' 

/joZcfe /or a// (x, z). 
Proof: 

By Ito's formula, we have: 

dP t {T) = exp (- J* R(Z s )ds^J (-R(Z t )P t (T)dt + dP t {T)) 

dP t {T) = ~9i<t>idt + Qi^dZt + 9 ~id\Z\t 

Setting the coefficient of dt to be 0 and we get the desired identity. □ 
Before moving on, we are now going to define some notation. 

i / \ , .dSAz) b 2 (z) d 2 <±>i(z) 
Mz) .= a(z)^->- + y 2* - R(z)Mz) 

Let Fk be the family of polynomials of degree at most k. That is: 

It follows immediately that F^ is closed under addition and subtraction and F& C i^+i- 

Theorem 2.2. Suppose the functions gi,<fri are given, if further suppose the functions gi(-) are 
linearly independent, we have that A,i{z) are uniquely determined and can be expressed as linear 
combinations of (j>i(z). Moreover in the polynomial interest rate model when (j>i(z) = z % , we have 
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Ai(z) eF n fori = 0,1,..., n. 
Proof: 

The identity of Thm 2.1 may be rewrite as: 

n n 

(*) = £ft0«04(s) 

i=0 i=0 

Fix any z, There are n + 1 unknowns, namely Ai(z). We can take n + 1 sample points Xj, j = 
0,1, ... ,n, such that we can rewrite (*) in the matrix form: 

Mx = c 

Where x = {A 0 (z), A n (z)) T , c = {Ya= 0 9i{xo)4>i{z), Ya=o 9i{x n )4>i{z)) T and the matrix has 
component Mij = gj{xi). Further by independency of fuctions gi, we may take M to be invertiblc 
and hence we have a unique solution x = M~ x c. Notice that the components of c are linear 
combinations of 4>i(z), so does the solution M _1 c. Finally since z was fixed arbitrarily, we get the 
required result. □ 

2.3. Necessary Conditions. 

Theorem 2.3 (Necessity). In polynomial models, if the coefficient functions gi{-) are linearly 
independent for i = 0, 1, . . . ,n then we have a(z), b 2 (z), R(z) are all polynomials and R{z) G F2, 
a{z) G F 3 , b 2 (z) G F 4 . Further if a(z) = £- = o a ^ R ( z ) = £i=0^ = £i=0^' We also 

have the following constraints on the coefficients: 

n(n — 1) , „ 
na 3 + — -b 4 - R 2 = 0 

(n-l)a 3+ (n - 1) 2 (n - 2) 6 4 -fl2 = 0 

n(n — 1) , „ 
na 2 + V ; 6 3 - i?i = 0 

Proo/: 

Recall the definition of and by Thm 2.2 we have Ai(z) G F n . Therefore: 

A 0 (z) = -R{z) 

Ai(z) = a{z) - zR{z) 

A 2 {z) = 2za{z) + b 2 (z) - z 2 R(z) 
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Hence the functions a(z), b 2 (z), R{z) are all polynomials in z. 
On the other hand, we have: 

A n {z) = nz n - l a{z) + ^ ~ ^ z n ~ 2 b 2 {z) - z n R(z) G F n 

Therefore we have: 

(1) nza(z) + nz n " x a{z) + ^ ~ % 2 {z) - z 2 R{z) G F 2 
Calculating A n _\ and A n _2 similarly we can get: 

(2) (n - l)za(z) + ^ - ~ 2) b 2 (z) - z 2 R(z) G F 3 

(3) (n - 2)za(z) + ( " ~ 2)( " ~ 3) fe 2 (z) - z 2 i?(z) £ F 4 
Since is closed under addition and substraction, (1) — (2) gives 

(4) za(z) + (n-l)b 2 (z) G F 3 
(2) - (3) gives 

(5) za(z) + (n- 2)b 2 (z) G F 4 

(4) - (5) gives 6 2 (z) G F 4 and hence za(z), z 2 R(z) G F 4 . i.e. a(z) G F 3 and R(z) G F 2 . 
For the relation between the coefficients of functions a, b 2 , R, we reconsider 

A n (z) = nz n ~ l a{z) + ^ ~ l \ n ' 2 b 2 {z) - z n R{z) 
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= nz 



1 E^M + ^f^- 2 ( j>< ) - * B ( 



\i=0 / \i=0 / \i=0 / 

Since A n (z) G F„, the coefficient of term z n+2 must vanishes and hence we have: 

n(n — 1) , „ 
na 3 + — — ^ — L h - R 2 = 0 

Similarly by setting the coefficient of z n+1 terms of A n (z), A n _i(z) to be 0, we get the remaining 
identity on coefficients. □ 



2.4. Sufficient Conditions. The above theorem provide necessary conditions on the functions 
a(z), b 2 (z), R(z) when searching for arbitrage-free polynomial interest rate models. On the other 
hand, these conditions are also part of sufficient conditions. 
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Theorem 2.4 (Sufficiency). Suppose a{z) = Y^=o a i zl > R( z ) = Yli=o ^ iZ% > b 2 (z) = Yli=Q^i zt are 
polynomials that satisfy the following conditions for some nature number n: 

n(n — 1) , „ 
nas + — ~b 4 - R 2 = 0 

(n-l)a 3+ (n - 1) 2 (n - 2) 6 4 -fl 2 = 0 
na 2 + n( ^ h -R 1 = 0 

And the factor process Z t with dZ t = a{Z t )dt + b(Z t )dWt is non- explosive. Then there exists a 
unique model Pt(T) = YA=o9i(T ~ O^t satisfying condition (A) where the spot interest rates is 
given by r t = R(Z t ). Moreover, if n < 3, the coefficients functions gi can be calculated explicitly 
and we have a factor interest rate model with explicit bond price. 

If in addition the spot interest rate process rt = R{Zt) is non-negative, then the discounted ZCB 

process Pt (T) is true martingale if and only if the factor process Zt is bounded. 

Proof: 



By theorem 2.1, if there exists functions gi(x) such that the corresponding identity holds with this 
particular choice of a(z), b 2 (z), R(z). Then we can form a polynomial model that satisfies condition 
(A). Therefore it remains to show that, we can always get a set of unique solution functions gi(x). 

The identity of Theorem 2.1 reads: 

i=0 i=0 V i=0 ' i=0 / 

Comparing the coefficients of z n+2 gives: 

/ n(n — 1) 
0 = g n I na 3 H b 4 - R 2 

Which is checked by the conditions of this theorem. Similarly comparing z n+1 terms gives: 

( n(n-l), \ /. . (n-l)(n-2), 
0 = g n [na 2 + — % - Ri) + g n -i (n - l)o 3 + 77 " fe 4 - #2 
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Which is also checked by the conditions of the theorem. In general, collecting coefficients of z l 
terms for i = 0, 1, . . . , n gives: 

(» + 2)(i + l), /. , i(i + l) ,\ ( 

9i = 9i+2~ 77 % + ft+i (i + l)ao + 0 % + # *ai + 0 - 62 - Ro 



(*) + ft-i ((i - l)a 2 + ( * ^ 2) b 3 - + 5,-2 ((i - 2)a 3 + ( * ^ 3) b 4 - i?2 

The formula (*) holds true provided we interpret terms g_2 , 9-i, 9n+i, 9n+2 to be 0. 

Having the above results, searching for an explicit polynomial interest rate model is equivalent 
to solve the coupled ODE 

(**) G(x) = SG(x) 

Where G(x) = (go(x), . . . , g n (x)) T and the (n + 1) x (n + 1) matrix S is given by: 

(» + 2)(i + l) 

*5j,j+2 - ^ 0 

S it i+i = (i + l)ao + V 2 61 

<Si,i = iai H - b 2 - Ro 

= (i - l)a 2 + (i ~ i y~ 2) b 3 - i?i 

S M _ 2 = (i - 2)a 3 + ( '~ 2) 2 ( '" 3) fe4 - R2 

with boundary condition G(0) = (1,0,..., 0) . We call the matrix S the information matrix 
because it provides all we need to solve the coupled ODE (**) and hence the identity in Theorem 
2.1. 

It remains to show that for any information matrix S, the coupled ODE (**) has a unique 
solution and we can form a unique polynomial model satisfies condition (A). This is true since we 
always have a unique solution to 1st order ODE with initial conditions. To be more specific, I will 
show that this is indeed the case for n = 1, i.e. the information matrix S is 2 x 2. For higher 
dimensions, the proof is similar. 

Case 1: suppose S has 2 distinct eigenvalue (maybe complex) Aj, with corresponding eigen- 
vector Uj. Then we have the eigenvectors must be linearly independent and the general solu- 
tion is then G(x) = ciw 1 e Ai:E + C2V 2 e X2X for some constants c\,C2- Applying boundary condition 
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G(0) = (1,0,..., 0) T , we have c\v x + C2U 2 = (1, 0) T . We can then solve for c\,C2 and this has a 
unique solution because vectors Vi,v 2 are linearly independent. 

Case 2: S has 2 repeated eigenvalue A with corresponding eigenvector v 1 . In this case, we can 
always find a vector v 2 that is independent of w 1 and that 

(S - \I)v 2 = 

where I is the identity matrix. In this case, the general solution to (**) is then 

G(x) = c x v x e Xx + c 2 (v 2 e Xx + v x xe Xx ) 

We now apply the boundary condition and still can get a unique solution. 

For Higher dimensions, the proof goes similarly. Especially, when n < 3, the information matrix 
S is at most 4x4. We know that formula exists for finding roots of quartic equation. Therefore, 
we can solve (**) explicitly and get a explicit bond price model. 

Notice that in general we can write the solution as 

G(x) = e Sx G(0) 

For the rest of this theorem, the discounted ZCB price process Pt(T) given by 

P t {T) = exp (- £ R(Z s )ds S j X>( T " l ) Z l 

is bounded if the factor process Zt is bounded. And we know that bounded local martingales are 
true martingales. 

On the other hand, if Pt(T) are true martingales for any T, then we have 
P t (T)=E[P T (T)\3? t ] 



9i(T - t)Z\ = E exp (- ^ R(Z s )ds) \& t 

i— n V \ J t J 



< 1 



Hence Zt must be bounded. □ 



Remark 2.3. Theorem 2.4 implies that we can freely choose the factor process and spot rate process 
as long as the conditions in Theorem 2.4 are satisfied. Then we are guaranteed to get a consistent 
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bond price model. Further, for quadratic and cubic models, we can always have an explicit bond 
price expression. 

Remark 2.4. Notice that the situation here is completely different as in the exponential polynomial 
case, where Filipovic [3] proves that the maximum degree is 2. Here we can choose the polynomial to 
an arbitrary large degree so long as the sufficient conditions stated in Theorem 2.4 holds. Actually 
this is the crucial reason that we can extend the model to the stochastic volatility case as we will 
see in section 3. 



2.5. Examples of Parametric Families. Therefore we can start by choosing some stochastic 
processes that capture some features of the spot rate and set up the corresponding bond price 
model depending on spot rate process. Here below are some possible choice for the factor process: 
(n = 2) 



Example 3 (Family 1) 



dZ t = aZ t (Z t - k)dt + y/pz t {k - Z t )dWt 
r t = 2aZ t 

with parameter a,f3,k> 0. The information matrix S of this family takes the form: 

( o n n \ 



0 



0 



V 



—2a —ka kj3 
0 -a -2ka-f3 J 

Remark 2.5. This choice of Z t provides non negative spot rate lives in the interval [0, 2ka\. Also 
we can see that the information matrix S have an eigenvalue 0, hence we can solve the remaining 
eigenvalues by solving a quadratic equation easily. However, the process might be absorbed at the 
boundary in finite time. 



Example 4 (Family 2). 

dZ t = a(p - Z t )dt + ^Z^dWt 
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with parameter a, (3 > 0. The information matrix S of this family takes the form: 



I 



\ 



0 a/3 0 
-1 -a 2a/3 
0 -1 -2a 



Remark 2.6. This family provides consistent models with factor process being spot rate process 
itself. Also the spot rate will be non-negative and fluctuating around some mean spot rate level (3. 
The parameter a measures the strength that pull the spot rate to the mean spot rate (3. There is no 
absorption issue in this family but the spot rate is not bounded above. 

Example 5 (Family 3). 



dZ t = a((3 - Z t )dt + y/Z t (k-Z t )(l-Z t )dW t 

n = z t 

with parameter a > 0 and 0 < f3 < k < I . The information matrix S of this family takes the form: 

( 0 a/3 0 "\ 

-1 -a 2af3 + kl 
0 -1 -2a -k-l 



S = 



V 



) 



Remark 2.7. This is a modification of family 2 that provides consistent models with factor process 
Zt being spot rate process r t itself. Also the spot rate will be in the interval [0, k] and fluctuating 
around some mean spot rate level (3. Moreover since the factor process is bounded, the discounted 
bond price provided by this family will be true martingales. 
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3. POLYNOMIAL STOCHASTIC VOLATILITY MODELS 



3.1. Model Setup. In this section, we would like to extend the models from last section and apply 
them to the stochastic volatility case. For simplicity in this section, we will set the spot interest 
rate r t = 0. i.e. there is no discounting in time. This assumption greatly simplify the notation 
and calculation without affecting the generality of our results provided the spot interest rates are 
deterministic. We will occasionally mention some changes on the results in the stochastic spot rate 
case where r t = R(Z t ). 

We model the power option price process Pt(T,8) as 



for some deterministic functions ki and ^. Again we assume that the factor process Z t follows: 



Similar to the interest rate model discussed above, we assume that the volatility process at also 
depends on the factor Zt, i.e. 



for some deterministic function h(z). 

Remark 3.1. If we set (fii(z) = z l , we recover the polynomial type model for the stochastic volatility 
case. 

The boundary condition Pt{T,0) = can be made to hold true by setting ko(0, 6) = 1, 
ki(0, 6) = 0 i = 1, 2, . . . , n and 4>o(z) = 1. By applying a similar argument as Remark 2.2, there 
is no loss of generality by assuming this type of boundary conditions on ki. Hence here follows, we 
will always use this type of boundary condition. 

3.2. No Arbitrage Condition. In the following subsection, we are going to establish the inter- 
connection between functions ki,4>i,a,b. And prove that under some constraints, Black-Scholes 
model is the only example. 




dZt = a{Z t )dt + b{Z t ).dWt 



We also assume that under pricing measure, the stock price process St follows: 



dS t = Stat.dWt 



a t = h(Z t ) 
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For consistency, we need: 

(B) : The processes P t (T, 6) are local martingale under some EMM Q V(T, 0) 

Remark 3.2. In the stochastic spot rate case, condition (B) will changed to require the processes 
Pt(T,8) are local martingales, where 

rt 



Pt(T, 6) = (exp - J R{Z s )d^j P t (T, 6) 



For ease of notation, we will introduce the following functions: 

A{x,e,z) :=f;fc^(*,#)<Mz) 

i=0 



B(x,0,z) :=J2ki(x,6) 



i=0 
n 



d<f>i(z) 
dz 



C(x,e,z) :=E^^)^P 

i=0 



D{x,0,z):=Y,0h{x,0) 

i=0 

For no arbitrage, we have: 

Theorem 3.1. Condition (B) holds if and only if 

(0 

n 

h{x, e)4>i(z) = \\h(z)\\ 2 A(x, 6, z) + a{z)B(x, 9, z) + \\b{z)\\ 2 C(x, 6, z) + b{z).h{z)D{x, 6, z)*) 

i=0 

holds for all (x, 6, z) 
Proof: 



Applying Ito's formula, we get: 



dSf = Sf [eh(Z t ).dW t + *0(0 - 1) \\h{z)\\ 2 dt 

d (j2 HT - t, e)UZt)) = jz (-kifr + ah d -p- + h t \\b\\ 2 0) dt + kA.dW t 

\i=0 ) i=0 ^ Z 2 ' Z 

Then we could use the product rule for stochastic process to get: 

dP t = sfd ( f>(r - t, e)UZt)\ + ( X>Cr - *, o)MZt)) ds e t + d ( J>(t - 1, e^zM ds 6 t 

\i=0 J \i=0 / \i=0 J 
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Finally we set the coefficient of dt to zero and get the desired identity (*). 
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Remark 3.3. In the stochastic spot rate case, we have an extra —R(z) Ya=o ki(x, 6)<j>i(z) term on 
the RHS of identity (*) 

Now if the functions h(x, 9) and <pi(z) are known, we would like to invert Theorem 3.1 and solve 
for functions a(z), b(z), h(z). This is handled by the following: 

Theorem 3.2. If for arbitrary z , the functions A(-,-,z), B(-,-,z), C(-,-,z) and D(-,-, z) are linearly 
independent, then the functions a{z), ||6(z)|| 2 , ||/i(z)|| 2 and b{z).h{z) are uniquely determined by 
functions ki(x,9) and (fii(z). 
Proof: 

Fix any z arbitrarily. There are 4 unknowns, namely a(z), ||6(z)|| 2 , ||/i(z)|| 2 and b{z).h{z) in the 
identity derived in Theorem 3.1. Also by linear independency, we can find 4 sample points (xj,0j) 
with % = 1, 2, 3, 4 such that the 4 by 4 matrix with ith row formed by vector 

(A(xi, 6i, z),B(xi, 9i, z), C(xi, 0,, z), D(x h 9 h z)) 

is invertible. Hence we can solve for a(z), \\b(z)\\ 2 , \\h(z)\\ 2 and b(z).h(z) and the solution is unique. 
Since z is fixed arbitrarily, we get the result. □ 

Remark 3.4. In fact as we are going to prove later, in the special case that 4>i{z) = z l , i.e. 
in polynomial model, we can say more about the functions a(z), \\b{z)\\ 2 , \\h(z)\\ 2 and b(z).h(z). 
Actually they are all polynomials in z. 

Now we define the following functions: 

Oi(z) := \\h(z)\\ 2 Mz) 

Theorem 3.3. if the functions e ( e ~ 1 ^ ki(x, 9), ki(x,9), 9ki(x,9) are linearly independent for i = 
0, ...,n. Then the functions cti(z) , Pi(z) , Ji(z) are uniquely determined. Moreover, they are linear 
combinations of (f>o{z) , ■ . ■ , <j> n (z) . 
Proof: 
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We can rewrite the identity established in theorem 3.1 as: 



n n 



0(0-1) 



2 



ki{x,9) + (3i{z)ki(x,6) + ~ /i (z)6k i (x,0) 



i=0 i=0 



By a similar argument as Theorem 2.2. We have that the unknowns are functions oti,/3i,^i and 
there are N = 3(n + 1) of them. By linear independency, we can take N sample points (xj,9j) and 
form an invertible matrix with its rows formed by functions ^ ki(x, 0), ki(x,9), 9ki(x,9) that 
evaluate at points (xj,9j). Therefore the functions cti(z), 0i(z), 7i(z) are uniquely determined. 

Moreover, notice that the matrix we formed above is a constant matrix, especially it is indepen- 
dent of z. Hence, when multiplied by its inverse matrix, we get the solution cti,f3i,~fi are purely 
linear combinations of 4>o{z), . ■ ■ , 4> n (z). □ 

For polynomial model, we then have the following corollary by setting 4>i{ z ) = z% '- 



Corollary 3.4. For polynomial models, when the independency holds, we have the functions 
ai(z),Pi(z),ii(z) e F n . 



3.3. Black-Scholes Model is the Only Example. We are now ready to prove the rather sur- 
prising result. 



Theorem 3.5. With the notation above, in polynomial models with factor process Z t , i.e. 



if the volatility process a t = h(Z t ) is a function of the factor process and we have the functions 
ki(x, 9), ki(x, 9), 9ki(x, 9) are linearly independent for i = 0, . . . , n. Then there is no consis- 
tent models except the trivial Black-Scholes model. 
Proof: 



n 



P t (T, 9) = Si ^ki(T- t, 9)Z\ 



i=0 



dS t = Sto-fdWt 
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By corollary 3.4, we have the functions ati(z), Pi(z), ji(z) are polynomials of degree < n. Espe- 
cially we have: 

ao(z) = \\h(z)\\ 2 hence ||/i(z)|| 2 is polynomial 
Pi(z) = a(z) polynomial 
faiz) = 2za(z) + ||6(z)|| 2 hence ||6(z)|| 2 is polynomial 
71(2) = b(z).h(z) polynomial 

At the same time, we have: 

a n {z) = z n \\h{z)\\ 2 a polynomial of degree < n 

But the right hand side cannot have degree < n unless the degree of ||/i(2i)|| 2 is 0, i.e. a constant, 
say a 2 . We can then define the following stochastic process: 

dW t = -a t .dW t 
a 

Notice that process Wt is a local martingale with quadratic variation 

d[W t ] = ^ htfdt = dt 

By Levy's characterisation of Brownian Motion, we conclude that Wt is also a Brownian Motion 
and we have dSt = StcrdWt- Hence with re-parameterizing the factor process Z± in terms of the 
new background Brownian motion Wt we have the resulting process St is then geometric Brownian 
Motion and we recover the Black-Scholes model. □ 

3.4. An Extension To Polynomial Model. In last section, we have seen that under some 
linearly independency constraints, Black-Scholes model is the only no-arbitrage model that makes 
power options price process Pt(T,9) have the form 

n 

P t (T, 6) = Sf ^h(T -t,6)Zl for all 6 

i=0 

Where 



dS t = S t h{Z t ).dW t 
dZt = a(Z t )dt + b(Zt).dW t 
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However, if we first fix 9, and then try to find a no-arbitrage polynomial model. We will be in the 
similar situation as interest rate case and hence we can have infinitely many examples. To be more 
specific, we will have infinitely many consistent models of the form 

n{6) 

P t (T,9) = S e t ^2k t (T-t,6)Zl(e) 

i=0 

dS t = S t h(Z t {9)).dW t 
dZ t {9) = a(Z t (0),9)dt + b(Z t (8), d).dW t 

That is for each specific value of 9, we choose a specific n{6) and a corresponding factor process 
Z t (9). 

So we may wonder can we tweak the original polynomial model such that we allow more no- 
arbitrage models instead of just Black-Scholes but don't need to choose a brand new factor process 
Zt every time we change the value of 6. The answer is yes and we propose a family 

D N = |^K = l,2,...,iV-l| 

When 6 G Djy, we set up our model (M) as 

n{6) 

P t (T, 9) = S 6 t HT - t, 9)Zl(D N ) 

i=0 

dS t = S t h(Z t {D N ),D N ).dW t 
dZ t (D N ) = a(Z t (D N ),D N )dt + b(Z t (D N ), D N ).dW t 

i.e. Now the factor process Zt is family dependent instead of 6 dependent. Hence we can use the 
same factor process Zt as long as we are in family D/v- For the following, in order to simplify the 
notation, we will drop Djy dependence on functions a,b,h and still write as a(z) , h(z) , b(z) . But 
bear in mind that we may choose them to depends on family Dn instead of 9. 

Remark 3.5. If we get a consistent model (M), we then get the exact modelling of price Pt(T,9) 
at 6 £ Dn- And hence we are able to approximate the remaining 9 G (0, 1). 

By applying Theorem 3.1, the no-arbitrage condition for our model proposed above is 

n(9) n{6) 

k(x, 6)<t>i{z) = k(x, 9)Bi(z, 9) 

i=0 i=0 
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Where 



= z HM^OII <Mz) + g?(z,#)^— + - ||6(z)|| 



2 11 v ; " rn y v ' ' dz 2 11 v ; " 8z 2 
<f>i(z) = z { 

d{z,9) = a(z) + 6b(z).h(z) 
By using a similar argument as Theorem 3.3, we have the following 

Theorem 3.6. If for fixed 9, the functions 0) are linearly independent for i = 0, 1, . . . , n. Then 
Bi(-,9) G F„(0) 

Remark 3.6. In the stochastic spot rate case, replace e ^~^ \\h{z)\\ 2 by 9( - e ~ 1 ^ \\h(z)\\ 2 — R(z) in 
the definition of Bi(z,9). 

3.5. Necessary Conditions and Sufficient Conditions. Similar to the case of interest rate, 
we get 

Theorem 3.7 (Necessity). In polynomial stochastic volatility model (M), if for any 9, the coefficient 
functions ki(-, 9) are linearly independent for i = 0, 1, . . . ,n then we have d(z, 9), \\h(z)\\ 2 , ||6(2:)|| 2 
are all polynomials in z and \\h(z)\\ 2 £ F 2 , d(z,9) G F 3 , \\b(z)\\ 2 £ F 4 . Further if d(z,0) = 
J2i=odi{0)z l , \\h(z)\\ = Yli=o hiZ % , \\b{z)\\ =Yli=Q^i z% - We also have the following constraints on 
the coefficients: 

Proof- 
Recall the definition of F>i and by Thm 3.6 we have Bi(z, 9) G F n ^ for each 9. Therefore: 
B 0 (z,9) = d ^ Y ^\\h(z)\\ 2 
B 1 (z,9) = d(z,9) + e ^l z \\h(z)\\ 2 
B 2 (z, 9) = 2zd(z, 9) + \\b{z)\\ 2 + fcil z 2 mz)ll 2 
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Hence the functions d(z,9), \\h(z)\\ 2 , ||6(z)|| 2 are all polynomials in z. 
On the other hand, we have: 

B n(e) (z,9) = n(9)z^d(z,9) + ^"f-V "-' \\Hz)f + z m e ^L \\h(z)f G F n(e) 
Therefore we have: 

(6) n(9)zd(z, 9) + mm^l Mz) f + , 2 fcD || M , )f e F2 
Calculating B n ^_ 1 and B n ^_ 2 similarly we can get: 

(7, (um - m ) + W-y-" ,«,),> + , ^ 3 

(8) („(«) - 2)^,0) + ( "«" - 2) 2 ( " ( " - 3) ||f>(.)|| 2 + z' 6 -^> \\ H z)f € F 4 
Since is closed under addition and substraction, (6) — (7) gives 

(9) zd(z,9) + (n(9)-l)\\b(z)\\ 2 GF 3 

(7) - (8) gives 

(10) zd(z,9) + (n(9)-2)\\b(z)\\ 2 G F 4 

(9) - (10) gives \\b(z)\\ 2 G F 4 and hence zd(z,ff),^^ \\h(z)\\ 2 G F 4 . i.e. d(z,5) G F 3 and 
||^)|| 2 GF 2 . 

For the relation between the coefficients of functions d(z, 8), \\h(z)\\ 2 , ||fr(2:)|| 2 , we reconsider 
B n(e) (z,9) = n(8)z^d(z,9) + ||6(,)|| 2 + ^)^_12 || M ,)f 

\i=0 / \i=0 / " Vi=0 

Since B n ^(z,9) G F n ^, the coefficient of term z n W+ 2 must vanishes and hence we have: 

„ W3W + "W("W-i) t4 + g^I) ft2 = o 

Similarly by setting the coefficient of z n ( e ) +1 terms of B n ^{z,9),B n ^_i{z,9) to be 0, we get the 
remaining identity on coefficients. □ 
Likewise with the same argument as theorem 2.4, we have 
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Theorem 3.8 (Sufficiency). Suppose d(z,9) = ELo^W*' l , \\H z )f = Ei=o ^ > \\ b ( z )\\ 2 = 
Yli=obiZ l are polynomials in z that satisfy the following conditions for some nature number n(9): 

,^ , ,m n(0)(n(0) - 1), 5(^-1), 
n{9)d 2 (9) + ^ ^ ^ 3 + V 2 % = 0 

^4n(i i/ie factor process Z t with dZ t = a(Z t )dt + b(Z t ).dWt is non-explosive. Then there exists a 
unique model Pt(T,9) = Sf Y^i=o ki{T — t)Zl satisfying condition (B) where the volatility process is 
given by at = h(Z t ). Moreover, in the case n{9) < 3, the coefficients functions ki can be calculated 
explicitly. 



3.6. Worked Examples. We will discuss some specific examples of model (M) in this subsection. 
We have already seen in section 3.3 that if ||/i(z)|| 2 is a constant, .i.e. polynomial with degree 
0. Then we will recover the Black-Scholes model. Hence in order to find interesting consistent 
examples, we must have ||/i(z)|| 2 at least a linear function, say ||/i(z)|| 2 = h\z + ho. By Theorem 

3.7, this corresponds to the case hi = 0 and hence we must have ds(9) = 0 and b^ = 0. Therefore 
by Theorem 3.8, the sufficient conditions becomes the factor process Z t is non-explosive and 

,m , ,„x n(9)(n(9) - 1) , 9(9 - 1) , 
(11) n{9)d 2 {9) + V A y >-b 3 + K - , h 1 = 0 

Given that 9 £ Djy, we must find suitable functions a, b, h that depend on N instead of depend on 
each single 9 and nature number n{9) that satisfy (11). 



Example 6 (Family 4). 



\\h{z)\\ 2 = 2N 2 z + hv 
\\b(z)\\ 2 = b 2 z 2 + b 1 z + bo 
b(z).h(z) = 0 

a(z) = z 2 + d\z + do 
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Where bi,di,ha are arbitrary constant that produce non-negative and non-explosive factor process 

z t . 

Proof: 

Since b(z).h(z) = 0, we have d(z,6) = a(z) = z 2 + d±z + do and hence function d depends on 
family Djy only with d 2 = 0. We can simply check that (11) holds true and by sufficient theorem 
3.8. No-arbitrage model is guaranteed. □ 

Remark 3.7. Notice that in above example, the functions a(z),b(z) are actually independent of 
family Djy. Hence the advantage is that we can use the same factor process Z t for any precision N. 
However the drawback is thatn{9) can be as large as which means the corresponding information 
matrix S can be large, i.e. ^ x Therefore solving the eigenvalues and eigenvectors of S is 
not easy. 

Remark 3.8. For the purpose of the non-negative and non-explosive factor process Z t , we can 
choose from the following parametric family 

dZ t = (Z t - ai ){Z t - a 2 )dt + y/ pz t (i - Zt)dW} 

Where the parameters are ft > 0 and 0 < ct\ < 7 < 02- This parametric family will produce factor 
process Z t G [0, 7] without any absorbing at boundary issue. 

Example 7 (Family 5). 

\\h(z)\\ 2 = cN 2 z + h 0 
\\b(z)\\ 2 = -cz 3 + b 2 z 2 + hz + b 0 
b(z).h(z) = 0 

a(z) = — - — cz 2 + d\z + do 

n( — ) = minji, N — i} 

Where c > 0 and bi, di, ho are arbitrary constant that produce non-negative and non-explosive factor 

process Zt- 

Proof: 



Simply check condition (11) and apply sufficiency Theorem. 



□ 



POLYNOMIAL MODELS FOR INTEREST RATES AND STOCHASTIC VOLATILITY 25 

Remark 3.9. Since the linear coefficient of \\h(z)\\ 2 is cN 2 , for any prescribed precision N, we 
could still choose c > 0 freely so that the norm of volatility process \\o- t \\ 2 = \\h(Z t )\\ 2 can be arbitrary 
linear function with positive linear term. Moreover we see that n{6) is at most y ; the complexity 
of the information matrix S is reduced significantly, i.e. at most y x y. 

Remark 3.10. For the purpose of the non-negative and non-explosive factor process Z t , we can 
choose from the following parametric family 

dZ t = ^Y^c(Z t - ai )(Z t - a 2 )dt + y/cZ t (Z t + /3)( 7 - Z t )dW} 

Where the parameters are (3 > 0 and 0 < ct\ < 7 < 02- This parametric family will produce factor 
process Z t G [0, 7] without any absorbing at boundary issue. 
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4. FURTHER ANALYSIS AND EXAMPLES 



In this section, we work through some examples that we purposed earlier with simulated graphic 
results. 

4.1. Stationary Distribution. To find the stationary distribution for a diffusion process Zt such 
that 

dZ t = a(Z t )dt + b(Z t )dW t 
We write p(t, x, y) be the transition density of process Z t in time t and let 

/oo 
g(x)p(t,x,y)dx 
-oo 



be the density at time t when the initial r$ have density g(x). Then we know that / satisfies the 
forward equation (Fokker-Planck equation). 

dt=-d-y {fa) + 2W ifb) 
If the initial density g(x) is stationary, then all following density will stay the same, i.e. 

f(t,y)=g(y) Vt 

Hence the stationary distribution g must satisfy the equation 

Solve the above equation gives 



(12) g{y) oc -L ( eX p J | 

4.2. Example on Interest Rate Case. Here we would like to have a close look at the Family 2 
that we purposed in section 2.5. For a recap, we construct a family of quadratic bond price models 
that depends on the spot rate process rt and takes the form: 

P t (T) =go(T-t) + 3i (T - t)r t + g 2 (T - t)rj 

Where the spot rate process r t comes from the family 

dr t = a(fi - r t )dt + JrfdWt 
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for some positive parameter a, (3. 

First notice that this family has an stationary distribution that can be calculated explicitly. By 
formula (12) with a(y) = a((3 — y) and b(y) = \/y^, we get that the stationary distribution R have 
a density function 



9r{t) oc -^exp 



2a a/3 



r > 0 



Which doesn't look very nice currently. A change of variable X = 1 / R reveals the density funtion 
of X is 



fx(x) = 9r(1/x) 



dr 
dx 



1 



oc xexp ( — a(3(x — — ) 2 



x>0 



To find the normalising constant of density fx, we integrate over [0, oo) and set the result to 1 
gives 

1 



fxix) = Cxexp ( -af3(x - 



x > 0 



C 



1 



2a/3 



e P +2 




Where $ denotes the cumulative standard normal distribution function. Hence the CDF of the 
stationary distribution is given explicitly by 

F[R <r]= P[X > -] 



I Cxexp ( — a/3(x — — ) 2 jdx 

Jl/r \ P J 



c 



1 f Rt 1 1 

exp {— ap{- 



2 } + 2 A /^ (v ^( 



1 1 

r ~ 0 



)) 



r > 0 



2a(3 "V (5' 

To solve the model completely, we need the information matrix S, in our case with above pa- 
rameter, we have 

( 0 1/20 0 ^ 
-1 -1 1/10 
0 -1 -2 



V 
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Solve the coefficient functions and the corresponding yield curve is given by 

Yt{T) _ .asm 

Here below, we provide simulated sample path of spot rate process r< and yield curve Yq(T) of 
Family 2, with parameter a = 1 and /3 = 5% . Recall that f3 can be interpret as the mean rate and 
we can see the path fluctuate between 0.05 and stays positive as predicted by Remark 2.6. 
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4.3. Sample Path on Stochastic Volatility Case. Here we simulate sample path of Example 
7 in section 3.6. We choose our parameters to be N = 100, c = 10 -4 , ho = 0, a\ = 0.02, a.\ = 0.06, 
P = 0 and 7 = 0.05. 
Hence our model reads 

dZ t = " X 2 10 (Z t - 0.02) (Z t - 0.06)dt + 0.01 x y/z?{0.05- Z t )dW? 
dS t = Sty/ZtdW? 

Where W l ,W 2 are independent Brownian motions. 

Remark 4.1. By Remark 3.7, we know that factor process Z t will live inside [0,0.05], which tells 
us that the volatility process at = h(Z t ) = \fZ~ t is around 20%. 

Here below, we provide a graph showing the sample path of stock price St with initial Sq = 1 
and a graph showing the implied volatility with maturity date T = 5, 10, 20 and strike K ranges 
from 0.5 to 2.5. 
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sample path of stock price S 




0 2 4 6 8 10 12 14 16 

Time t 
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5. Appendix A 

Here we will define the forward variance model and prove that under some linear independency 
conditions, if the forward variance process is polynomial function of factor process Z%, then the 
maximum degree is 2. 

We work under the pricing EMM Q and model the stock price process by dSt = Stcr t .dWt. We 
parameterize the process {Pt(T, 9)} by 

P t (T,0) = ^ e -^(i-e)/ t T /(t, s; e)d s 
Where df(t, T; 6) = B(t, T; 6).dW t + A(t, T; 9)dt 



Theorem 5.1 (No arbitrage conditions). If we have 

(13) A(t, T; 6) = B(t, T; 6).{ 6{l ~ 6) £ B(t, s; 6)ds - 6a t ) 
and the process f(t, T; 9) satisfies 

(14) f(t,t;9) = \\a t \\ 2 V(t,0) 

Then the price process of power options P t (T,9) are local martingales for all (T,9). 
Proof: 



Let M t = 9logS t - \9{l - 9) jf f(t, s; 9). Then P t (T, 9) = e Mt is a local martingale if and only 
if the drift term of dM t + ^d[M] 4 vanishes, where [M] denotes the quadratic variation of M. By 
applying Ito's formula, we have: 

dM t = 9jdS t - ^d[S] t + l -9(l - 9) f{t, t; 9)dt - l -9{\ - 9) jT df(t, s; 9)ds 

= (9a t - h(l - 9) jT Bit, s; 9)ds).dW t - (^0(1 - 6) jT A(t, s; 9)ds + )f \\a t \\ 2 )dt 



d[M] t 



(^9 2 \\a t \\ 2 -9 2 (l-9)a t . B(t,s;9)ds + h 2 (l-9) 2 B(t,s;9)ds 



dt 
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Hence the dt term of dM t + \d[M} t is 



-0(1-0) / A(t,s;e)ds- 1 e 2 (l-e)a t . f B(t,s;6)ds+ 1 e 2 (l-e) 2 [ B{t,s;0)d. 

2 Jt 2 J t 8 7t 

= --9 2 {l-9) 2 f B(t,s;9). [ B{t,u-6)duds + \e 2 {l-6) 2 [ B(t,s;6)d 

4 J t Jt ° Jt 



which vanishes by symmetry. □ 

Remark 5.1. Changing into Musiela [6] notation, i.e. Setting f t (x,9) = f(t,t + x;9) we can 
recover the model we introduced in Section 1. 

For general factor models, if we require the forward variance to have the form 

f t {x,0) = f(t,t + x;9)=g(x,9,Z t ) 

Hence 

f(t,T;9) = g(T-t,9,Z t ) 

Where the factor process Z t follows 

dZ t = a(Z t )dt + b(Z t ).dWt 

for some continuous deterministic functions a(z), b(z),g(x, 0, z). Moreover, if we insists the process 
at depends also on the factor process Z t , say 

a t = h(Z t ) 

Then by applying Theorem 5.1, we can get the no arbitrage conditions for general factor model. 

Theorem 5.2 (No arbitrage conditions for general factor models). Suppose the functions a(z), 
b(z), g(x,9,z), h(z) satisfy the following 

9x = \ INI 2 (gzz ~ 0(1 - 0)g z jf g z (y, 0, z)dz^j + (9b.h + a)g z 
g(0,9,z) = \\h(z)\\ 2 V(M) 
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Then we have an arbitrage- free factor model for forward variance by setting: 

f(t,T;9)=g(T-t,9,Z t ) 
a t = h(Z t ) 
dZ t = a(Z t )dt + b(Z t ).dW t 

Proof: 



By applying Theorem 5.1, we have 

g(0,9,Z t ) = f(t,t;9) = \\a t \\ 2 = \\h(Z t )f 
Applying Ito's formula on f(t,T; 9) gives: 

df(t, T; 9) = -g x dt + g z dZ t + -g zz d[Z] t 

(15) = bg z .dW t + (ag z + - ||6|| 2 g zz - g x )dt evaluated at point (T - t, d, Z t ) 

Hence by reading off the coefficient dWt term, we must have B(t, T; 6) = b(Z t )g z (T — t, 9, Z t ) and 
by (13), we have: 

A(t, T; 9) = g z b. ( * heta ^ Z 9 ) j T 9z{s _ tj 9j Z t )b(Z t )ds - Oh{Z t )j 

(16) = \\b\\ 2 9(1 ~ 6) g z £ 1 g z (x,9,Z t )dx-8g z h.b 

Equating (15) and (16) gives the required results. □ 



Integrate the no arbitrage condition of Theorem 5.2 with respect to the first argument from 0 
to x. Writing 



G(x,0,z) = [ g(y,9,z)dy 
Jo 



We have: 



(17) G x = \ \\bf (G zz - e ^ Y e ±G 2 z ) + (9b.h + a)g z + \\h(z) 

(18) G{0,9,z) = 0 Vz 
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For polynomial forward variance models, G(x,9,z) = X^=o^( x '^) z * ^ or some deterministic 
functions k{. Then (17) takes the form: 

Theorem 5.3. We can generate an arbitrage- free polynomial forward variance model if 

t "t^y = \ www 2 (x>(M)|£ - 6 ^(±H X ,e)%A 

i=0 \i=0 i=0 / 



(19) +(eb(z).h(z) + a(z))(J2ki(x,e)^) + \\h(z^ 2 

i=0 

holds for all (x, 0, z) and ki(0, 9) = 0 for all 9 and i 



Define 



1 2 dz* dzi 

aij(z) = aji (z) := - ||&(*)|| 

n/i(z) := b(z).h(z) — 



We can rewrite (19) as 



n r\, n n 

J2 ^{x, e)z i = J2 Hx, 9)Pi{z) + OH*, oH(z) 

i=0 i=0 i=0 

n 

(20) - ^(1 - °) k i(x, °)W x i 0)<*ij{z) + IIM 

i,j=0 



z)\\ 2 



Theorem 5.4. // the functions l,ki(x,9),9ki(x,9),9(l — 9)ki(x,9)kj(x,9), 0 < % < j < n, are 

1 1 2 

linearly independent. Then aij,{3i,ji, \\h\\ are uniquely determined and are polynomials of z with 

degree < n 

Proof: 



Set C = 1 + (n + 1) + (n + 1) + — — — - be the total number of unknown functions ctij(z), fii(z), 
ji(z), \\h(z)\\ 2 . Similar to the proof of Theorem 2.2, we can take C sample point (x m , 9 m ) for m = 
1,2, ... ,C and form an invertible CxC matrix whose m-th row consists of 1 , ki (x m , 9 m ) , 9 m hi (x m , 9 m ) , 
9 m (l — 9 m )ki(x m , 9 m )kj(x m , 9 m ). Hence the solution is unique. 
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In addition, we have the matrix we formed above is independent of z. Therefore its inverse is 
also independent of z and the solutions must be linear combination of z°, . . . , z n . □ 

We now prove our main result on polynomial forward variance model which is similar as in the 
interest model case [3]: 

Theorem 5.5 (Maximal Degree in Polynomial Forward Variance Model). Given any polynomial 
model, if the functions l,ki(x,9),9ki(x,9),9(l — 9)ki(x,9)kj(x,9), 1 < i < j < n, are linearly 
independent. Then there is no arbitrage- free model for n > 3 
Proof: 

Recall the definition of a, we have: 

/x 1,,, ,„2^9z 

a ^ = ~J b ^ dzdz 
= \\\b(z)\\ 2 

t ^ 1 ml/ M.2 dz n dz n 

(*) =\\\b(z)\\ 2 n 2 z 2n - 2 

But by Theorem 5.4 we have both a nn {z) and ||6(2;)|| 2 are polynomial with degree at most n. But 
the right hand side of (*) cannot have degree at most n unless n < 2. □ 
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